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ABSTRACT 

Aims. UltraCompact Hn regions are signposts of high mass star formation. Since high-mass star formation occurs in clusters, one 
expects to find even earlier phases of massive star formation in the vicinity of UltraCompact Hn regions. Here, we study the amount 
of deuteration and depletion toward pre/protocluster clumps found in a wide-field (10 x 10 arcmin) census of clouds in 32 massive 
star-forming regions that are known to harbour UCHn regions. 

Methods. We determine the column density of NH 3 , NH 2 D, CO, H 13 CN, and HC 15 N lines. We used the (J,K) = (1, 1) and (2,2) 
inversion transitions of NH 3 to constrain the gas temperatures. 

Results. We find that 65% of the observed sources have strong NH2D emission and more than 50% of the sources exhibit a high 
degree of deuteration, (0.1 < NH2D/NH3 < 0.7), 0.7 being the highest observed deuteration of NH3 reported to date. Our search for 
NHD 2 in two sources did not result in a detection. The enhancement in deuteration coincides with moderate CO depletion onto dust 
grains. There is no evidence of a correlation between the two processes, though an underlying correlation cannot be ruled out as the 
depletion factor is very likely to be only a lower limit. Based on simultaneously observed H 13 CN and HC 15 N (7=1-0) lines, we derive 
a high abundance ratio of H 13 CN to HC l3 N, which might indicate anomalous ratios of C and N isotopes relative to those derived 
toward the local ISM. 

Conclusions. We find CO depletion and high deuteration towards cold cores in massive star forming regions. Therefore, these are 
good candidates for sources at the early phases of massive star formation. While our sensitive upper limits on NHD 2 do not prove 
the predictions of the gas-phase and grain chemistry models wrong, an enhancement of ~ 10 4 over the cosmic D/H ratio from NH 2 D 
warrants explanation. 

Key words. ISM: molecules -chemistry - deuteration - Stars: formation 

1. Introduction earlier phases of massive star formation, and the raw material 

out of whi ch massive stars or sta r clusters form, in the vicinity 

The earliest phases of high mass star formation are still poorly of UCHns (Thompson et al. 2006h. 

understood. The prevalent theory is that massive stars are born in ^ has p rompted us t0 embark upon a p r0 g ram t0 sea rch 

dense clumps within giant molecular clouds (GMCs), where the for massive p re /p ro tocluster clumps by taking a wide-field (10 x 

presence of Ultracompact HII regions (UCHII) identifies mas- 10 arcmin) cen sus of clouds in 32 massive star-forming re - 

sive stars which have already begun to ionize their surround- gior ^ harbouring UCHn reg i ons dWood & Churchwelll lT989). 

ings. Observations dCesaroni et al. || 1994l | 1998|) have shown that Qur program is known as SCAMPS (the SCUB A Massive 

these are often associated with warm (T > 100 K), compact(< p re /p ro tocluster core Survey: [Thomson et al.ll2005L Thompson 

0.1 pc) and very dense («h, > 10 cm ) cores known as hot et a j 2007 in prep ) 

molecular cores (HMCs). An evolutionary sequence has been ~ We discovered a multitude f secondary, cold condensations 

suggested based on recent observations ( Codella et al.| |2004l and wamer clumps that might contain heavily embedded mas- 

| Beuther& Shepherd ! | 2005|) with the HMC stage immediately siye protostars Many of the sources are seen as dark patches 

preceding the formation of UCHn regions. But the stage be- in MSX i nfrared im ages of the region, infrared dark clouds 

fore hot cores - the precluster and early protocluster phases of (TRDCs; Menten et al.l200lh . As a result these clumps must have 

massive star formation - have not yet been studied. Since many dust temp eratures below 30 K (as evidenced by MIR upper lim- 

UCHn regions are located in clusters, one expects to find even ks) and haye masses of a fgw lQQ tQ 1000 Mq Thefje dumps 

thus could be in a colder pre-protostellar phase. 

Send offprint requests to: tpillai@cfa.harvard.edu In order to study the physical and chemical conditions of 

* Present address: Harvard-Smithsonian Centre for Astrophysics these clumps, we carried out a multi-wavelength survey toward 
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them in various molecular tracers. In this study, we also in- 
cluded 9 sources which were selected on the basis of the MIR ex- 
tinction (IRDCs) and followed up later with submillimeter dust 
continuum emissio n and millimeter rotational lines of H2CO 
dCarev et al.lll998l) . We have recently rep orted a study of am- 
monia in this sample (iPillai et al.l 120061) . Here, we report on 
our study using ammonia to probe the temperature and deuter- 
ation in the clumps and CO to estimate depletion of neutral 
molecules from the gas phase. Depletion has proved to be a 
goo d tracer of low-m ass pre-stellar cores (Bac mann et alj|2003l 
and lCrapsi et al.ll2005l) . and NH3 is an important tool in measur- 
ing the physical conditions in molecular clouds and can be used 
to confirm the low temperatu res and high densities required for 
depletion dHo & Towne si [19831) . Since only the lowest NH3 en- 
ergy levels are expected to be populated for cold dark clouds 
{T < 20 K), their physical conditions can be probed using the 
(1,1) and (2,2) inversion transitions in the metastable (J, K) rota- 
tional levels of ammonia. Radiative transitions between different 
Zf-ladders are forbidden, therefore the lowest levels are popu- 
lated only via collisions. The optical depth can be determined 
from the ratio of the hyperfine satellites. Thus, the population of 
the different levels can be estimated and hence the temperature 
of the gas determined. 

NH3 is observed to be an excellent tracer of the dense 
gas where many ot her molecules would have heavily depleted 
dTafalla et alj 120041). In addition, recent chemical models re- 
veal that NH3 (and also N2H + ), does not deplete from the gas 
phase for the densities observed in dark clouds (< 10 6 cm -3 ; 
Ber gin & Langerill997l) . However in the light of recent experi- 
ments of the desorption of CO and N2 at temperatures relevant to 
dense cores, the relative difference between the derived CO and 
N2 binding energies is found to be signific antly less than tha t 
currently adopted in astrochemical models (Ob erg et al.l 12005). 
Therefore, the observed non-depletion of NH3 and N2H + during 
the early formation of a massive cluster remains a "mystery". 

High densities in the earliest phase are expec ted to enhance 
deple tion of molecules (mainly CO) onto grains (Bacmann et al. 
120031) . The freeze out of abundant gas phase molecules 
along with the low temperatures (< 20 K) preva lent in these 
clump s leads to a high degree of de uteration dFlower et alj 
2005LlRoueffet aU l2005l lHatchellll200l IShah & WoottenT hoO 1 , 
Saito et al. 2000). Hence, the estimation of deuteration and de- 
pletion in clumps can be used as a chemical and physical 
chronometer. 

In this paper we investigate the NH3 temperature and deuter- 
ation in our candidate pre/protocluster cores using observations 
with the Effelsberg 100m and IRAM 30m telescopes. In Sect. 
|2l we describe our single dish observations. We also discuss 
the data reduction and present the derived physical parame- 
ters for different molecules. In Sect. [3j we show that candidate 
pre/protocluster clumps exhibit very high deuteration and deple- 
tion. The fractionation of NH3 is compared with depletion of CO 
from the gas phase. The possible origin of deuterated species 
and the variation of deuteration with gas temperature is also dis- 
cussed. HCN was also observed as part of the 30m observations 
and we conclude the discussion with our interpretation of the 
HCN isotopic abundance with galactic distance. 

2. Observations & Results 

In this section we present the results of single dish observations 
of NH 3 , NH 2 D, CO and HCN lines on a sample of 23 SCAMPS 
sources and 9 IRDCs. The brightest compact submm sources 
with no MIR counterpart were selected for these observations 



from our original catalogue of the SCUBA fields. The parame- 
ters of the molecular transitions covered, along with the observ- 
ing setups are given in Table Q] In Table |2j we give the source 
list. The derived line parameters and the column densities are 
presented in Tables 3-8. 

2.1. IRAM 30m observations 

The NH 2 D, C ls O, H 13 CN and HC 15 N observations were made 
at the IRAM 30m telescopeQ on Pico Veleta in August 2003 to- 
ward the 32 candidate pre/protocluster clumps given in Table [2] 
We performed 9 point maps with 10" spacing around the peak 
to be able to compare, independently of the beam, lines at differ- 
ent frequencies. All line parameters listed are averaged over the 
map. 

Individual positions were observed with an integration time 
of 2 minutes per position in the position switch mode. The 
receivers A100, B100, A230 and B230 were tuned to 86.13, 
109.975, 219.56 and 231.32 GHz respectively. 

The VESPA autocorrelator was used at the backend, with a 
spectral resolution of 40 kHz and 320 kHz at 100 and 230 GHz 
respectively. We had average summer weather conditions with 
a maximum system temperature of 248 K at 110 GHz. Towards 
selected sources we also have C n O (J = 2-1) and NHD2 (lio- 
loi) observations. The half power beam width (HPBW) of the 
30m telescope is 22" at 1 10 GHz and 1 1" at 230 GHz. The main 
beam efficiency at these frequencies is 0.75 and 0.52, respec- 
tively. 

2.2. Effelsberg 1 00m observations 

We observed the NH3 (1,1) and (2,2) transitions with the 
Effelsberg 100m telescope in April 2004 for the 23 sources 
listed in Table|2] Two of the sources were re-observed because of 
baseline problems with the receiver. With the AK 8192 backend, 
we were able to observe the (1,1), (2,2), (3,3) and (4,4) tran- 
sitions in both polarisations simultaneously using the K-band 
receiver. With 8 subunits of 10 MHz bandwidth, the resulting 
spectral resolution was as 0.25 kms -1 . The beam at the NH3 
frequencies was 40". The observations were performed in the 
frequency switching mode. Pointing was checked at hourly in- 
tervals by continuum scans on G10.62. We estimate the pointing 
to be accurate to within 6 ". The pointing scans were used for 
the absolute flux calibra tion. The NH3 o bservations toward the 9 
IRDCs are reported in lPillai et al1(l2006l) . Note that for G12.19- 
0.12 and G29.97-0.05 we report the observations carried out in 
October 2002 with a spectral resolution of w 1 km s _1 . 

2.3. Results of the NH3 and NH2D single dish observations 

The ratio of the brightness temperatures of the NH3 (1,1) 
and (2,2) transitions, along with the optical depth, can be 
used to estimate the rotational temperature. For temperatures 
<20 K, which are typical for cold dark clouds, the rota- 
tional temperature closely follow s the gas kinetic tempera- 
ture ( Walmslev & Ungerechts 1983). Like NH3, its isotopologue 
NH2D also has hyperfines which allows the estimation of optical 
depth and hence the column density. Thus the fractionation ratio 

1 IRAM is an international institute for research in millimeter astron- 
omy. IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). 

2 Based on observations with the 100m telescope of the MPIfR (Max- 
Planck-Institut fr Radioastronomie) at Effelsberg). 
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can be estimated assuming that NH2D and NH3 are co-spatial 
and hence have the same gas temperature. 

Out of the 32 sources that were observed, NH2D with hy- 
perfines were detected in 22 sources with S/N ratio > 5cr while 
we have a 100% NH 3 (1,1) and (2,2) detection. While we de- 
tected the (3,3) lines in most of the sources, the (4,4) lines were 
not detected. The fits to the hyperfines for NH3 and NH2D were 
done using CLASS methods NH3 (1,1) and HFS respectively 
dForveille et al.l 1989b . The line parameters from the resultant fits 
are given in Table IA.2I and Table IA.3I For those sources with a 
high uncertainty in the NH2D optical depth, the main line is fit- 
ted with a single Gaussian and the integrated intensity is quoted. 
The spectra are shown in Fig. IA.ll 

The basic physical parameters, namely the rotational tem- 
perature, the kinetic temperature and ammonia column density 
were derived using t he standard formulation for NH3 spectra 
(Bach iller et aflfl987h . The expressions used to estimate the col- 
umn densities from the radiative transfer equations for all other 
molecules including NH2D are given in Appendix A. The de- 
duced physical parameters are given in Table [3] The uncertain- 
ties given in brackets are formal errors obtained by a Gaussian 
error propagation. We have derived the NH2D column density 
solely based on the ortho transition. We assume that the ortho 
and para transitions are in LTE. 

For the few cases with significant detection of the para state 
we find that ratio of the integrated line intensity ( J TmbcIv) is 
approximately 3 consistent with a thermalised ortho to para ratio 
of 3. 

For the column density determination of NH2D, we have av- 
eraged over the 9 points of the map to compare to the NH3 re- 
sults. Essentially the resolution of the NH2D observations has 
been degraded to that of the NH3 observations. A caveat is that 
a map obtained with 10" spacing at the IRAM 30m does not 
completely sample the 40" beam of the Effelsberg telescope. 
However, it is very unlikely that there is significant emission 
outside the beam of the Effelsberg telescope (i.e. its FWHM) 
because a)we find that the brightest peak is always at the center 
of the map, b) box averaging delivers a similar result to a 1 D 
Gaussian weighted average c) the peak dust emission positions 
have been chosen as the centre coordinates for these observa- 
tions. 

The rotational temperatures are in a range from 11 to 23 K, 
with ammonia column densities from 1 to 5xl0 15 cirT 3 . We ob- 
tain [NH2D/NH3] ratios from 0.004 - 0.66. A deuteration factor 
of 66% in N H 3 is much large r than the largest reported value so 
far,~33% bv lHatchelll (12003). In 10 out of the 22 sources with 
NH2D detection, we obtain abundance ratios <0.02 while the 
rest of the sources have a high degree of deuteration (>13%). 
These 10 sources with low ratios show optically thin NH2D 
emission. 

In Fig. Q] the correlation between the velocities and line 
widths of the NH2D and NH3 emission is shown. Velocities and 
line widths are clearly different in some sources, with the NH2D 
line widths being smaller in most cases. This could be either due 
to NH2D and NH3 tracing different regions or the slightly differ- 
ent beams at the two frequencies, in which case, clumping might 
contribute to a larger line width. It is interesting to note that the 
critical density of the reported NH2D transition is a factor 47 
higher than that of the corresponding NH3 transition, assuming 
similar collisional rates. However, a more likely explanation is 
that NH3 could be originating from the more disturbed proto- 
cluster environment, while NH2D might be tracing the gas in 
precluster cores (Pillai et al. in prep.). 



2.4. Results of the CO observations 

The C ls O (1 — > 0) and (2 — > 1) observations were performed 
simultaneously with NFbD (1 1 1 - lo,i)- Figure [A2\ shows the 
spectra toward all sources. Note that there are secondary fea- 
tures in the C ls O spectra for a few sources, likely to be line of 
sight components not seen in dense gas like NH3. In Table I A. 41 
we list the observed line parameters based on Gaussian fits, and, 
in Table [4] the column density and excitation temperature esti- 
mates (see Appendix A) and the C 18 abundances. The C 18 
excitation temperature has been derived from the C 18 J = 1-0 
and 7 = 2-1 line ratios. 

Based on the dust continuum and CO observations for the 
observed sample, the depletion can be studied for the first time 
on a sample of massive clumps. Recent observational studies on 
the condensation of CO in low mass prestellar cores reveal that 
in the initial cold and dense evolutionary stage, CO is heavily 
deple ted onto the dust grains dBacmann et al.ll2002l ISavva et all 
2003). The depletion factor 77 is defined as 



Ni 



c 18 o 



[Nh 2 ] 
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Nc 18 o 
[Nh 2 ] 
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The C 18 column density, Nc<no> is derived from the observed 
integrated intensity (see Appendix A) and assumimg that the gas 
temperature is equal to the excitation temperature derived from 

the C ls O line ratios. 

The effective H2 column density (Launhardt 1996) is calcu- 
lated from 



JV(H 2 ) 



6.2 x 10 16 S V A 3 e 



zo 



(2) 



where S Y is the flux density in Jy/beam. A the wavelength in 
fim, Z/Zo is the metalicity relative to the solar metalicity(we 
assume 1) and 9 is the FWHM of the beam at wavelength 
A. We used an opacity K m of 1.85 cm 2 /g of dust at 850 fim 
for dust grains with thick ice mantles and gas density n(H) = 
10 6 cm 3 (Ossenkop f~& Henningll9 94). The column density was 
derived after smoothing the 850 //m dust continuum emission to 
the 20" resolution of the C ls O data (Thompson et al. 2007 in 
prep). 

In these dense clumps it is possible that the C 18 line may 
be optically thick and hence derived column densities could be a 
lower limit. In order to investigate this possibility, we observed 
the rarer isotopomer C I7 toward two of our sources, G35.19- 
1.73 and G8 1.74+0. 59 as shown in Fig. [2] The line parameters 
are given in Table|5] In Sect. l3.1.Tl we discuss how C n O is used 
as a more reliable tracer of the optical depth. 

2.5. Results of the H 13 CN and HC 15 N observations 

We also observed the isotopomeric species of HCN, H 13 CN 
and HC I5 N in our 3 mm setup with the IRAM 30m telescope. 
H I3 CN was detected in almost all sources except G19.30P1/P2 
and G10. 15-0.34. In Fig. IA.3I we show the spectra toward all 
sources where both H 13 CN and HC 15 N are detected with suffi- 
cient S/N (> 3 cr). The H I3 CN line has hyperfines, which allows 
to estimate the optical depth. However, we find that H 13 CN is 
optically thin in almost all cases. In two sources with moderate 
optical depths, the uncertainty in t is too large to derive mean- 
ingful constraints on the column density. The line parameters are 
given in Table [6] along with their column densities. The column 
density is computed using the gas temperature derived from NH3 
as the excitation temperature (see Appendix A). 
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Fig. 1. Left panel: Comparison of the LSR velocities of NH3 and NH2D. The median spectral resolution after Hanning smoothing 
is marked in the upper left corner of the panel. Right panel: Comparison of NH3 and NFbD line widths. Points with no y error bars 
are those for which the line width has been fixed based on the hfs fit to the optically thin emission, to derive the NH2D integrated 
intensity. 
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Fig.2. C 17 (J = 2-1) spectra toward G35. 19-1.73 and 
G8 1.74+0.59 observed with the IRAM 30m telescope. 

3. Analysis and Discussion 

3.1. CO Depletion 

The derived C ls O abundances are a f actor of 2 -12 lowe r 
than the canonica l value of 1.7-10 -7 dFrerking et al.l 119821) . 
IWu & Yangl d2005l) find similar depletion factors for the infrared 
dark cloud G79.2+0.38, another of the clumps in the cloud in the 
present study. 



the ISM. As discussed in lKramer et al.l d!999l) . the measured ra- 
tios of the integrated intensities of the C ls O and C n O lines can 
then be used to determine the C 18 optical depth. 

Due to the limited observing time, we observed the C 17 



18,17 



line only in the two brightest C O sources. If and r 
denotes the ratio of the integrated intensities of the C 18 and 
C n O lines and the optical depth of C 18 as derived from the 
ratio, then 



R 



18,17 



1 



•exp(-r! 8 ' 17 ) 



l-exp(-r; 8 - 17 /A(18,17)) 



(3) 



The assumptions involved in deriving Eqf3] are that both iso- 
topomers have the same excitation temperatures and beam filling 
factors. 



For G35. 19-1.73, we find that r 



18,17 



1.35, which im- 



plies that the correction factor for the column density given by 
t/(1 — exp(-r)) = 1.8, is clearly smaller than the observed de- 
pletion. However, for G81.74+0.59, the ratio itself/? 1817 ~ 6, 
suggesting that Eq.Q is not valid anymore. It is unlikely to be 
due to different excitation conditions or the extent of emission of 
C ls O and C I7 0. It is indeed possible that there are real variations 
in the ratio of 18 to n O, a possibility which needs to be investi- 
gated. Additionally, the observed relatively small observed peak 
brightness temperatures for C 18 exclude high optical depth in 
this line. Therefore, to the extent that we assume G35.19-l.73 
to be representative of our sample, we may conclude that C ls O 
optically depth is "moderate". 



3.1.1. Is C ls O optically thin? 

As mentioned in Sect. l2.4l the CO column density as determined 
from C 18 transitions might be underestimated if C 18 is opti- 
cally thick. However, the rarer isotopomer C 17 is expected to 
be optically thin. The canonical value of the relative abundan ce 
of C 18 w.r.t C n O, A(18, 17), is 3.65 dWilson & Rood|[T994l) in 



3.2. NH 3 deuteration 

How do we account for the high deuteration we observe in 
these pre/protostellar clump s? There are two main pathways 
to br i ng about deuteration dRodgers & Charnlevl 1200 ll |Millar 
l2002l lMillarll2003l iRoueff et :al]|2005l) : a) gas-phase reactions, 
b) production of deuterium bearing molecules on grain surfaces. 
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Accretion of neutrals onto the dust grains that would other- 
wise destroy H2D + enhances deuteration. Up to now all these 
processes have been used to explain fractionation in different 
sources. 

In pre/protostellar cores, ion-molecule exchange reactions 
prevalent at low temperatures coupled with the depletion of CO 
from gas phase is preferred over grain produ ction of highly 
deuterated species at te mperatures of 20 K dShah & Woottenl 
l20QllCaselli et al.ll2003h . The freeze-out of molecules from the 
gas onto the dust grains, particularly that of heavy molecules 
like CO in the cold initial phase h as been predicted by chem- 
ical models dFlower et al.) 120051 iRodgers & CharnleyJ _[2001. 
Roberts & Millar! l2000allb[ iBrow n & Millar 1989bl Iwatsonl 
19761) . These predictions have been validated by observations of 
molec ular freeze out onto dust grains in low- mass starless dense 
cores dCaselli et al.| [l999. Tafalla et al. 2002). The primary frac- 
tionation reaction that dominates at low temperatures (< 20 K) 
is 

+ HD -> H 2 D + + H 2 + AEi (4) 

Neutral molecules like CO can destroy H 2 D + , thereby lowering 
the deuterium enhancement. Therefore, the depletion of CO, the 
second most abundant molecule, from the gas phase can lead to 
an enhancement in the [H 2 D + ] / [Hi] ratio and thereby the molec- 
ular D/H ratios. 

iRoberts & Millar! d2000bl) show that at 10 K, accretion of 
neutrals onto the dust grains, especially CO, leads to the forma- 
tion of doubly deuterated molecules such as NHD 2 and D 2 CO. 
Based on these arguments, we expect to see a correlation be- 
tween deuteration and CO depletion. In Fig. [3] we compare the 
NH3 fractionation and the degree of CO depletion (depletion fac- 
tor). As seen from the correlation plot, there is a large scatter in 
the values and we do not see any obvious trend of high depletion 
and deuteration. Nevertheless, the main result is that we find very 
high deuterium fractionation and CO depletion. 

Depletion could well occur on scales much smaller than 
the beam for sources at the distances typical for these 
sources(several kpc). In that case, the CO depletion that is mea- 
sured as an average over the telescope beam is an underesti- 
mate of the true depletion. The NH 2 D emission, on the other 
hand, is probably dominated by the dense clumps. This could 
explain why the CO depletion does not appear to track the 
deuteration in Fig. |3] The uncertainties on the CO depletion 
factors are quite large, mainly due to the uncertainty in the 
H 2 column density estimate from the 850/im dust continuum. 
The larg est uncertainty is a factor of 4 in the dust opacity at 
850/um (|O ssenkop f & Henning| [l994, Kriigel & Siebenmorgen 
119941 iDraine & Leelll984l) . But we note that with the choice of 
the Ossenkopf & Henning opaciti es, the derived column d ensi- 
ties are already lower than for e.g. Savage & Mathisl dl979 ) dust 
properties. The two studies other than ours tha t directly com- 
pare C O de pletion and d e uteri um fraction are bv lBacmann et al.l 
d2003l) and lCrapsietail (120051) . Bacmann et al. use the D 2 CO 
to H 2 CO ratios to determine the deuteration. They claim to find 
a correlation between D 2 CO to H 2 CO ratios and depletion and 
argue that the presence of a significant amount of O (which is 
also an important H^ destroyer) in the gas phase still undepleted 
might be responsible for the large scatter observed. Crapsi et al. 
observe N 2 H + and N 2 D + toward 3 1 low-mass starless cores and 
find a good correlation between deuterium fractionation and CO 
depletion. 

If /(X) is the fractional abundance of species X relative to 
H 2 and R(XD) is the abundance ratio of XD relative to XH, then 
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Fig. 3. Comparison of NH3 fractionation and CO depletion fac- 
tor for the dense clumps in our sample. The error bar on CO 
depletion factor is quite large mainly due to the uncertain dust 
opacities adopted (see Sect. I3. It . 

under the assumption of steady-st ate and equ ating formation and 
destruction in Eq.d4jl, one finds dMillaril2003l) 

R(U 2 D + )^SuMT)f(HD). (5) 

where S h,d+ (T) is a function of the different rate coefficients, 
(forward, backward and dissociative recombination rates) given 
by 

kif 

Sh 2 d+(T) = - — ■ — 7777:- (6) 

k ir + «i/(e) + ££ M /(M) 

Here, k\f, k\ r are the forward and reverse rate coefficients 
such that k\ r = kifexp(-AEi/T) and a\ oc 1/ a/T is the dis- 
sociative recombination rate, /(e) is the electron fraction and 
X^m is the rate coefficient of H 2 D + with species M. AE\ is the 
exothermal energy (~ 220 K), also given in Eq.©. 

At low temperatures (T < 20 K) observed in cold cores, the 
primary fractionation reaction is dominated by the forward reac- 
tion in Eq.dUi. At higher temperatures, the reverse reaction be- 
comes important in removing H 2 D + very rapidly from the gas 
phase. Consequently, the decrease in primary fractionation pro- 
duces a corresponding decrease in secondary fractionation re- 
actions responsible for deuterium enhancements in molecules. 
Therefore, one naively expects to find an anti-correlation be- 
tween temperature and deuterium fractionation. In Fig. |H we 
plot the NH3 fractionation against the temperature for the 
high mass candidate pre/protostellar clumps (this paper) and 
for pre/protostellar clumps in the lower mass regime from the 
literature. The late st gas-phase predictions for [NH 2 D]/[NH3] 
dRoueff et alj|20 05) are also plotted. This model assumes a den- 
sity of n(H 2 ) = 10 5 crrT 3 , carbon and oxygen depletion factors of 
5 and 15 respectively while Nitrogen is kept constant. Our data 
points nicely fill in the missing observed data points between 10 
- 20 K. Clearly the agreement with the model is moderate. 

The levels of deuteration do fall off rapidl y with increas- 
ing temperature above 25 K (Roueff et al. 2005), but at the low 
temperatures shown here the dependence is small. The Roueff 
et al. models employ fixed H 2 densities of around 10 5 cm -3 , 
which would be typical averages over the 40" beam for our 
sources. However the NH 2 D emission will be enhanced at the 
peaks, where higher densities and higher depletion would be ex- 
pected to produce higher levels of deuteration. Depletion, which 
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increases over time roughly with the increase in density, has a 
strong influence on deuteration. The final levels of depletion and 
deute ration depend on the accretion history (eg. iFlower et al.l 
120051) . but it is clear that higher d ensities or lon g er tim escales 
result in higher levels of both, eg. iRoberts et al. (2003) derive 
[NH 2 D]/[NH 3 ] ratios of 0.4-0.8 for densities of 3xl0 6 cm 4 . 

However, most of the low [NH2D]/[NH3] ratios and in- 
terestingly the high temperatures are from the sources where 
NH2D column densities were estimated from the integrated in- 
tensity of the line and therefore are really lower limits on the 
[NH2D]/[NH3] ratios. The exact temperature which we measure 
for a clump of course depends on the combination of hot/cold 
gas in the beam, which can explain the scatter between 15 - 18 K 
where in some cases we are still seeing a lot of deuteration and 
in others it has already diminished. 

On average, the fractionation observed in sources in our sam- 
ple is very high. However, note that few sources at tempera- 
tures above 15 K show an order of magnitude lower deutera- 
tion than the model predictions, while others show extremely 
high deuterations, much higher than those found in low mass 
pre/protostellar cores. 

There could be three main reasons for the very low fraction- 
ation estimated for a few sources. First, one of the main assump- 
tion in deriving the [NH3]/[NH2D] ratio is that the filling factors 
for both molecular transitions are the same. Any deviation from 
this assumption might result in a discrepant ratio. Second, these 
sources might be relatively more evolved. In such a case, the ki- 
netic temperatures derived from the NH3 (1,1) and (2,2) lines 
might only be a lower limit. In such cases where the gas tem- 
perature is roughly higher than 20 K, the rotational temperature 
does not represent the real gas temperature. Hence the column 
density (deuteration) determined may be under estimated. Third, 
the sources could be chemically young, hence the timescale to 
reac h the high deuteration in these sources is larger than their 
age. lTafalla & Santiago! (2004) recently discussed such a chem- 
ically young low mass core. It is also possible that NH2D might 
be tracing different regions from NH3 and that the temperature 
derived from NH3 may not be the temperature of NH2D. 

3.2.1 . Upper limits from NHD 2 

A simple model for grain surface formation of mu l tiply de uter- 
ated molecules is considered by iBrown & Millar! <[l989a) and 
predicts the abundances of deuterated species, in order to dif- 
ferentiate between the gas and grain chemistry. According to the 
Brown and Millar model prediction, the abundances of deuter- 
ated species scale as [NHD 2 ]/[NH 3 ]= ±([NH 2 D]/[NH 3 ]) 2 for 
grain surface for mation. 

Alternatively, Roberts & Millar (2000a) propose a gas-phase 
chemistry where the effects of the freeze out of gas phase species 
onto grains is included and find an enhancement in the fraction- 
ation of both singly and doubly deuterated species. 

Very deep observations of the NHD2 lio-loi line for two 
sources G35.19-l.73 and G8 1.74+0.59, did not result in any 
detection at an rms of w 20 mK. Based on the 3 <x upper limits, 
we then derive [NHD 2 ]/[NH 3 ] (see Appendix A) and find that 
([NHD 2 ]/[NH 3 ]) observed < (0.02) G35 .i9; (0.02) G81 . 74 
(i([NH 2 D]/[NH 3 ]) 2 ) observed < (0.06) G35 .i9;(0.04) G81 . v4 
([NHD 2 ]/[NH 2 D]) observed < (0.05) G35 .i9;(0.05)g8i.74 

The grain-chemistry model prediction is based purely on 
probability arguments assuming that the fractionation is propor- 
tional to gas-phase atomic D/H ratio and that there are three dif- 
ferent pathways for inserting the D atom into the chain of re- 
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Fig. 4. Plot of NH3 fractionation versus gas temperature as de- 
rived from NH 3 . The dashed li ne is the latest gas phase model 
predictions (Roueff et al. 200H). The filled and u nfilled squares 
mark the SCAMPS sources and iHatchelll (I2003I) respectively; 
filled and unfilled trian gles, the values found to war d low mass 
pre/pr otostellar cores by Shah & Woottenl d200ll) a ndlSaito et al. 
(2000) respectively and the unfilled circles mark iRoueff et al. 
(2005) sources. 



actions that form multiply deuterated NH3. In spite of the sim- 
plicity of the model, the prediction is close to the derived upper 
limits in both sources. 

How do the obtained upper limits compare with the predic- 
tions of recent gas-phase deuterium chemistry models? In Fig. [5] 
we compare the results of model calculation s with our observe d 
upper limits for the ratio of NHD 2 to NH 2 D dRoueff et al]2005l) . 
We have used their "model 2" which gives the highest molecu- 
lar fractional abundances for the N-bearing species. In this con- 
text it is also interesting to note that while Roueff et al. find 
that their model underpredicts the NHD2 to NH2D ratio in low 
mass prestellar cores (LDN 134N,LDN 1689N, Barnard 1, LDN 
1544), it agrees well with the observed ratio in our sample. 

Hence, we cannot reject active grain or gas-phase chem- 
istry on the basis of the observed NHD2 upper limits. While 
a gas-phase model along with condensation of neutral species 
onto grain mantles adequately explains moderate deuteration 
(NH2D/NH3) observed towards some sources, the very high 
fractionation we derive remains unexplained. Although we do 
not understand the grain-surface chemistry, particularly in this 
case involving deuterated ammonia, a naive extrapolation from 
the surface-chemistry involving H, D, C and O and a high grain- 
surface D/H ratio mi ght explain the sig nificant enhancement of 
deuterated ammonia (ICaselli et al.ll2004T) However, the efficacy 
of this model is strongly dependent on the gas-phase atomic D/H 
ratio, which in the extreme case of our sample implies an en- 
hancement of ss 10 4 over the cosmic D/H ratio. 

3.2.2. Comparison with deuteration in other star forming 
regions 

Although the exact evolutionary status of the sources (i.e. 
pre/protostellar) in our sample cannot be determined with cer- 
tainty, there are several observational evidences that substantiate 
their high mass pre/protostellar nature. The large line broadening 
(1 < Av < 5 kms~') observed in all molecular tracers, together 
with the high H2 column density (7Y(H2) > 10 23 crrT 2 ) is char- 
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Fig. 5. Comparison of obtained upper limits for G3 5.19-l.73 and 
G8 1.7 4+0. 59 with the gas phase model predictions (Roue ff et alj 
I2005L their model 2). 



distance from the galactic centre (kpc) 



Fig. 6. H 13 CN and HC 15 N abundance ratio as a function of the 
distance from the galactic centre. The results are consistent with 
a single value of Ri so = 17 shown as the dashed line. 



acteristic of massive star forming regions. These sources have 
masses between a few hundred to a few thousand solar masses. 
Their source-averaged H2 densities are a few 10 5 cm 4 - unlikely 
for low-mass objects. Moreover, the physical properties of these 
cores are similar to th ose associated with UCHn regions (e.g. 
iThompson et al.ll2006l) . Therefore we expect them to be forming 
(or on the verge of forming) high-mass stars. The exact nature of 
these sources will be discussed in a future paper (Hatchell et al. 
in prep.). 

The deuteration of a few tens of percent in ammonia is com- 
parable to the highest which have been measured in the inter- 
stellar medium, which occur in low-mass prestellar cores. The 
related ratio [N2D + l/r N2H + 1 exceeds 20 % in a few prestellar 
cores including L1544 dCrapsi et al.ll2004l) . High deuterated am- 
monia abundances have also been measured in early-stage pro- 
tostars: the rNH2Dl/[NH3| ratio reaches over 30% in protostars 
in Perseus (Hatchell 2003), and detections of tri ply deuterated 
ammonia (|va n der Tak et al. 2002; [Lis et al.120021) and methanol 
( Pari se et al.1 20041) have all been in the environments of proto- 
stars, though deuteration of ~ 5% is more c ommon ([Hatchell 
l2003HSaito et al.l2000tlShah & Woottenl2001l) . The high deuter- 
ation of our sample supports the idea that these are the high-mass 
equivalent of prestellar or protostellar cores. 

Hot molecular cores, believed to be the early protostellar 
stage of high mass star formation and potentially a later evo- 
lutionary stage of our sources, have a low molecular deuter- 
ation fractionation ~ 10~ 3 , consistent with ice formation at a 
higher temperature of 60 - 80 K and little gas-phase deuteration , 
as expected at hig h temperatures (Hatch ell et al.1 IT9981 [l999; 
iRoberts et ai1 l2002). High [NH 2 D]/[NH 3 ] ratios are observed in 
the O rion region b ut only in the compact ridge and not the hot 
core dTurnerill9 90). If our candidate precluster cores are precur- 
sors of hot cores, then the chemistry must radically alter during 
the later evolution of the hot core, perhaps with continuing ac- 
cretion at higher temperatures during the early protostellar phase 
dominating the ice production. 

3.3. Variations in HCN isotopic abundance 
Let us define the abundance ratio, ^; so , as 

Also - ^ u ^ ■ (7) 



j?i so is found to vary w i th the distance from the Galactic center 
(IWilson & Roodlll994l IWielen & Wilson! 1 19971) . The observed 
trend can be explain ed by Galactic chemical evolution. Other 
observational studie s (llkeda et al.l2002llLanger & Penziasll993l 
iDahmen et al.l l 1 995 ) have found in addition a source-to-source 
variation of the isotopic ratios among a group of clouds located 
at nearly the same distance from the galactic centre. 

The isotopomeric species of HCN, H 13 CN and HC 15 N have 
been often observed in their 7=1- transi t ion toward various 
local molecular c louds (Dahmen et al.ll 19951 iHirota et"ai1 1998. 
Iked a et aill2002l) to estimate the ratio ,Rj so . These lines can be 
observed with the same receiver and hence the ratio derived will 
be independent of calibration errors. Moreover, the optical depth 
can be measured from the three hyperfine components of the 
H 13 CN line. As mentioned in Sect. [53] H 13 CN is optically thin 
in almost all cases and hence the observed ratios must indicate 
the abundance ratio 7?; so , since the column density estimates are 
free from bias due to photon trapping effects. 

In Fig. [6] we plot the abundance ratio of H 13 CN and HC 15 N 
as a function of distance from the galactic centre for our sources. 
The ratio of the abundances does not seem to correlate with the 
distance. However, the formal error bars are too large to be able 
to make any strong argument against a correlation of /?; so with 
galactocentric distance. 

If we adopt the canonical values for the 12 C/ 13 C of 77 and 
14 N/ 15 N of 450 toward the local interstellar medium (ISM), 
then the implied H 13 C N to HC 15 N abundance ratio is ~ 6 
(IWilson & Roodl[l994[ their Table 4). From our data, we find 
that most of the values tend to be higher, consistent with a mean 
value of 17. This ratio is intermediate between that found toward 
the galactic centre (~ 30), an extreme high mass star-forming 
region and the local ISM (~ 6) where typically low mass sta r 
formation pervades (see again Table 4 of I Wilson & Roodll 19941) . 

The SCAMPS sources are in the close vicinity of massive 
star forming regions and hence the energeti c UV radiation coul d 
cause a isotope-selective photo-destruction ( Schi lke et al.l 1992b . 
Chemical fractionation is also likely to play an important role in 
the observed high values. The high deuterium fractionation we 
observed must be equally affecting the abundances of H I3 CN 
and HC 15 N. Therefore, if at all, chemical fractionation were to 
play a role, it must be attributed to an anomalous fractionation 
in the isotopes of either N or C. A high ratio implies either a 
low 15 N fractionation or I3 C fractionation in HCN. As noted 



8 



Pillai et al. 2007: Deuteration and Depletion in High Mass Clumps 



earlier in this sectio n, the carbon isotop ic ratio is known to vary 
across the Galaxy dWilson & Roodl 19941) . The ratio is influenced 
mainly by a) the self-shielding property of CO, b) 13 CO selec- 
tive isotopic fractionation, and c) stellar winds from stars more 
mass ive than the s un, which favours the production of 13 C over 
12 C (Rriigel 2003). The deviation from the canonical value in 
the local ISM could be therefore due to a complex chemistry or 
environment in these high-mass star forming regions. 

Although the number of assumptions involved is consider- 
able, in particular that the two species may not be exactly tracing 
the same volume of gas as evidenced by the differences in their 
line widths, the ratio we derive clearly indicates the differences 
in the properties of high mass and low mass star forming regions. 



4. Conclusion 

In this paper, we report the study of physical and chemical prop- 
erties of a new sample of high-mass pre/pro tocluster clumps. We 
have observed 32 sources in the NH 3 (1,1) and (2,2), NH2D (In— 
l m ), C ls O(7 = 1-0 and 2-1), H 13 CN and HC 15 N (J = 1-0) 
transitions. We have clear detections of NH2D in 22 sources and 
100% detection in NH 3 . 

Our observations suggest large deuterium enhancements (up 
to 66%), the largest reported so far. We also derived the amount 
of molecular depletion by comparing the C ls O column density 
with the H2 column density derived from dust continuum obser- 
vations, and we find that the degree of CO depletion is a factor 2 
- 12 higher compared to the canonical value for its abundance. 
The derived abundance ratio of H 13 CN to HC 15 N is indicative of 
a very high 13 C/ 12 C ratio, discrepant from that found toward the 
local ISM. These properties clearly reflect the complex chem- 
istry in regions of high mass star formation. 

This study shows that depletion and high deuteration exist 
towards massive cold cores in massive star forming regions and 
makes them promising candidates for the early phases of mas- 
sive star formation. In a subsequent paper, we will discuss the 
spatial distribution of NH2D in two sources observed with high 
angular resolution, confirming the high deuteration found from 
this study. 



Appendix A: Radiative Transfer Equations used for 
Column Density Determination 

We used the following expressions to determine the column den- 
sity of the different molecules. The column density is given by, 



In the optically thin case, 



3hso 



J(T ex )Q(T ex ) T Av , 



(A.l) 



where, Av is the linewidth, so is the dielectric permittivity, Sfi 2 
is the line strength multiplied by the dipole moment along the 
molecular g-axis, t is the optical depth, h is the Planck constant, 
r ex is the excitation temperature and Q(T ex ) is the partition func- 
tion. 

Here, 7(T ex ) is defined as 



/(Tex) = 



exp(E u /kT ex ) 
exp{hv/kT ex ) - 1 



(A.2) 



where E u is the upper energy level, k is the Boltzmann constant 
and v is the frequency of observed transition. We assume that 
the excitation temperature is the same for all rotational levels. 
Equation dA.lt is valid for the optically thick case. 



Ntat = 



3hep J(T ex )Q(T ex )W 
2n 2 Sfi 2 J v (T ex ) - J v (2.7) ' 



(A.3) 



where W is the integrated intensity and J V (T) is defined as 

1 (T -. _ hv/k 
Jv '' 1 > ~ exp(hv/kT a )-\ ' 

The NH2D partition function is determined by considering 
the contribution of the different energy levels from / = to J = 
2, while the metastable levels from (J,K- 1, 1) to (J,K = 3, 3) 
have been considered for NH3 . 

For other molecular species, the partition function Q(T ex ) at 
temperature T ex is estimated as Q(T ex ) = aT^, where a and /3 are 
the best fit parameters from a fit to the partition function obtained 
from JPL catalogue at different excitation temperatures from 10 
- 300 K. 

In Table lA. fl we give the dipole moments we used for differ- 
ent molecules and the partition function, Q(T ex ) at temperature 
r ex = 15K. 

By fitting the main and the hyperfine components of the (1,1) 
line and the main component of the (2,2) line, we obtain the 
rotational temperature. The kinetic temperature and NH3 col- 
umn density, hav e been derived using the standard formulation 
for NH 3 spectra dBachiller et al.lll987l) . The NH 3 column den- 
sity has been estimated assuming that the level populations are 
thermalised, i.e.,r ex = T rot . A similar assumption holds for our 
estimation of NH2D column densities. 

The upper limit on the NHD2 column density is estimated 
from the integrated intensity. This quantity was measured by in- 
tegrating over a velocity range determined from the NH2D line 
width around the LSR velocity of the source. The lcr uncer- 
tainty on this value is estimated from the rms of the noise in 
the channels where no line emission is expected, and, is given 
by VMv«ct. 
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Table 5. C n O (J = 2-1) line parameters 



Source 


J r M Bdv 


Av 




Kkms 1 


kms 


G35. 19-1.73 


7.1 (0.1) 


3.3 (0.1) 


G81.74+0.59 


5.5 (0.1) 


3.3 (0.1) 



Table A.l. Molecular parameters used to estimate N t 



Transition 


s& 


2(15) b 




D 




NH 2 D l„-loi (para) 


11.9 


47. 8 C 


NHD 2 lio-loi (ortho) 


0.722 


50.2 d 


H 13 CN J = 1-0, F = 2-1 


14.84 


22.8 e 


HC 15 N J = 1-0 


8.909 


7.6 e 


C 18 J = 1-0 


0.012 


6.1 e 


C ls O J = 2-1 


0.024 


6.1 c 



a S is the line strength and /i„ is the dipole moment along the molecular 
g-axis. 

b the partition function corresponding to a temperature of 15 K, the typ- 
ical temperature of the sources in the present sample. 
the partition function obtained by the summation over the contribution 
of the different energy levels from J = to J = 2. 
d the partition fu nction based on data from the CDMS catalog 
jMiiUer et alj200l) . 

e the p artition function based on data from the JPL catalog ( Pickett et al. 
119981) . 
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Table 1. Parameters of observed rotational transitions 



Species 


Transition 


A, (K) 


V (MrlZ ) 


JNH3 


(7, A = i, 1) 


22. /U 


23094.490 


iNri3 


{J, A — Z, L) 


AT CO 




NH 2 D (para) 




16.55 


85926.3 


NH 2 D (ortho) 


1 1 1 — lot 


15.98 


110153.6 


NHD 2 (ortho) 


lio-loi 


13.33 


110812.9 


NHD 2 (para) 


lio-loi 


13.09 


110896.7 


c 18 o 


J = 1-0 


0.0 


109782.1734 


c 18 o 


J = 2-1 


5.27 


219560.3568 


c 17 o 


J = 2-1 


5.39 


224714.3850 


H 13 CN 


J = l_0, F = 2-1 


0.0 


86340.184 


HC 15 N 


J = 1-0 


0.0 


86054.961 



Note: a E! is the lower energy level of the transition. 

Table 2. Positions and velocities of the observed sources. 1RDC positions are taken from Carey et al. 1998 



Source 


R.A.(2000) 


Dec.(2000) 


vlsr [km/s] 


SCAMPS 


G8.13+0.25 


18:02:55.69 


-21:47:46.7 


19.4 


G8.68-0.37 


18:06:23.24 


-21:37:14.1 


35.2/38.1 


G8.71-0.37 


18:06:26.51 


-21:35:46.6 


38.1 


G10.21-0.31 


18:09:20.63 


-20:15:04.5 


12.8 


G10.21-0.32 


18:09:24.52 


-20:15:41.4 


12.8 


G10. 15-0.34 


18:09:21.38 


-20:19:32.8 


12.8 


G10.61-0.33 


18:10:15.62 


-19:54:46.6 


74.0 


G12.19-0.12 


18:12:41.67 


-18:25:19.8 


27.6 


G13.18+0.06 


18:14:00.92 


-17:28:41.2 


51.6 


G15.01-0.67 


18:20:21.22 


-16:12:42.2 


26.2 


G15.03-0.65 


18:20:18.80 


-16:11:22.6 


26.2 


G15.01-0.69 


18:20:24.22 


-16:13:22.8 


26.2 


G18.17-0.30 


18:25:07.53 


-13:14:32.7 


54.9 


G18.21-0.34 


18:25:21.55 


-13:13:39.5 


54.9 


G23.41-0.23 


18:34:45.74 


-08:34:21.2 


104.2 


G23. 42-0.23 


18:34:48.16 


-08:33:56.1 


104.2 


G23.44-0.18 


18:34:39.25 


-08:31:36.2 


104.2 


G27.29+0.15 


18:40:34.70 


-04:57:18.1 


26.0 


G27.31+0.18 


18:40:32.45 


-04:55:03.8 


26.0 


G29.97-0.05 


18:46:12.25 


-02:39:05.9 


100 


G34.8 1-0.28 


18:53:20.148 


+01:28:31.5 


59.0 


G35. 19-1.73 


19:01:45.45 


+01:13:21.5 


42.4 


G81.74+0.59 


20:39:00.37 


+42:24:36.6 


-3.1 


IRDCs 


Gil. 11-0.12 PI 


18:10:34.04 


-19:21:4 


29.2 


Gil. 11-0.12 P3 


18:10:07.25 


-19:27:2 


29.2 


Gil. 11-0.12 P4 


18:10:07.25 


-19:28:4 


29.2 


G19.30+0.07 PI 


18:25:58.14 


-12:04:4 


26.3 


G19.30+0.07 P2 


18:25:52.69 


-12:04:4 


26.3 


G28.34+0.06 PI 


18:42:50.9 


-04:03:1 


78.4 


G28.34+0.06 P2 


18:42:52.4 


-03:59:5 


78.4 


G33.71-0.01 


18:52:53.81 


+00:41:06.4 


104.2 


G79.34+0.33 


20:32:26.20 


+40:19:40.9 


0.1 
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Table 3. NH2D and NH3 column densities, rotational temperature and fractionation 
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source 


1 rot 


Wnh 3 b 




[NH 2 D/NH 3 ] d 




K 


10 15 (10 14 ) cm" 2 


10" cm" 2 




G8. 13+0.25 


18.0 


(3.0) 


1.8 (0.4) 


1.7 (0.1 ) 


1.0 (0.2) 




17.4 


(2.3) 


1.6 (0.3) 


2.4 (0.1 ) 


1.5 (0.3 ) 


G8.68-0.37 


12.1 


(0.8) 


1.3 (0.1) 


39.4 (12.4) 


30.3 (9.8 ) 


G8.71-0.37 


13.3 


(1.0) 


1.5 (0.2) 


59.7 (10.6) 


39.8 (8.9) 


G10.21-0.31 


15.4 


(1.5) 


2.1 (0.3) 


70.9 (29.9) 


33.7(15.0) 


G10.21-0.32 


14.1 


(0.7) 


2.1 (0.1) 


51.3 (22.3) 


24.4 (10.7) 


G10.61-0.33 


17.8 


(1.5) 


2.3 (0.2) 


1.0 (0.1 ) 


0.4 (0.1 ) 


G11.11-0.12P1 


12.5 


(1.2) 


1.4 (0.2) 


54.0 (23.5) 


38.6(17.7) 


G11.11-0.12P4 


10.6 


(1.6) 


1.8(0.5) 


1.7 (0.3) 


0.9 (0.3 ) 


G12.19-0.12 


12.6 


(0.3) 


4.4 (0.2) 


5.6 (0.2 ) 


1.3 (0.1 ) 


G13. 18+0.06 


17.4 


(1.3) 


2.8 (0.2) 


88.9 (23.2) 


31.7 (8.6) 


G18.17-0.30 


15.8 


(1.2) 


1.8 (0.2) 


119.7 (31.5) 


66.5 (19.0) 


G18.21-0.34 


15.5 


(1.2) 


1.7 (0.2) 


47.3 (19.7) 


27.8 (12.0) 


G23.41-0.23 


18.8 


(1.8) 


2.5 (0.3) 


77.7 (51.3) 


31.1 (20.8) 


G23.44-0.18 


23.3 


(1.5) 


5.2 (0.4) 


5.0 (0.2 ) 


1.0 (0.1 ) 


G27.29+0.15 


17.8 


(1.7) 


3.0 (0.3) 


1.4 (0.2) 


0.5 (0.1 ) 


G28.34+0.06P2 


15.7 


(2.0) 


2.2 (0.3) 


4.1 (0.2) 


1.9(0.3) 


G29.97-0.05 


13.7 


(1.4) 


2.2 (0.3) 


67.7 (15.2) 


30.8 (8.1 ) 


G33.71-0.01 


19.5 


(4.6) 


3.0 (0.8) 


1.4 (0.1 ) 


0.5 (0.1 ) 


G34.8 1-0.28 


16.6 


(1.0) 


1.6 (0.1) 


2.7 (0.2 ) 


1.7 (0.2) 


G35.19-l.73 


17.0 


(0.9) 


1.8 (0.1) 


74.0(11.0) 


41.1 (6.5 ) 


G79.34+0.33 


13.7 


(1.0) 


1.4 (0.1) 


17.9(15.9) 


12.8(11.4) 


G8 1.74+0.59 


18.4 


(1.1) 


2.7 (0.2) 


92.9 (17.6) 


34.4 (7.0 ) 



11 Rotational temperature derived from NH 3 observations. 
b NH 3 column density. 
c NH 2 D column density. 

d Ratio of NH 3 and NH 2 D column densities in percentage. 

Note: Similar notations are used for other molecules elsewhere in this paper. 
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Table 4. C 18 column density and abundance from (J = 1-0) and (J = 2-1) transitions 



Source 




N c i& 


N(H 2 ) a 


g — 


rf 




K 


lO'^lO 15 ) cnr- 


10 2J cm" 2 


io- 8 




G8. 13+0.25 


11.8 


1.6 (0.1) 


2.6 


6.2 


2.7 


G8.68-0.37 


8.5 


2.4 (0.9) 


7.9 


3.0 


5.6 


G8.71-0.37 


4.2 


1.8 (0.5) 


3.5 


5.2 


3.2 


G10. 15-0.34 


19.1 


2.1 (0.4) 


3.8 


5.4 


3.1 


G10.21-0.31 


8.4 


1.4 (0.2) 


2.2 


6.4 


2.7 


G10.21-0.32 


8.0 


2.6 (0.3) 


7.4 


3.6 


4.8 


G10.61-0.33 


9.1 


1.0 (0.3) 


2.7 


3.8 


4.5 


G11.11-0.12P1 


6.4 


0.6 (0.1) 


1.1 


5.5 


3.1 


G11.11-0.12P3 


6.9 


0.5 (0.3) 


0.8 


6.2 


2.7 


G11.11-0.12P4 


6.6 


0.5 (0.2) 


1.4 


3.5 


4.8 


G12.19-0.12 


11.6 


1.2 (0.2) 


2.4 


5.1 


3.3 


G13. 18+0.06 


15.7 


1.2 (0.6) 


6.3 


1.8 


9.3 


G15.01-0.67 


32.5 


6.0 (0.4) 


16.6 


3.6 


4.7 


G15.01-0.69 


52.2 


4.5 (0.6) 


11.6 


3.9 


4.4 


G15.03-0.65 


51.7 


5.5 (0.6) 


14.2 


3.9 


4.4 


G18.17-0.30 


14.3 


1.5 (0.8) 


3.3 


4.6 


3.7 


G18.21-0.34 


10.9 


0.7 (0.4) 


2.5 


3.0 


5.8 


G23.41-0.23 


8.4 


3.4 (0.8) 


9.3 


3.6 


4.7 


G23.42-0.23 


9.0 


3.1 (0.3) 


3.8 


8.0 


2.1 


G23.44-0.18 


9.2 


3.0 (0.5) 


7.7 


3.9 


4.3 


G27.29+0.15 


8.8 


1.2 (0.2) 


1.7 


6.9 


2.5 


G27.31+0.18 


8.3 


0.5 (0.2) 


1.7 


2.7 


6.3 


G28.34+0.06P1 


4.8 


1.3 (0.5) 


1.5 


8.6 


2.0 


G28.34+0.06P2 


8.7 


1.4 (0.5) 


7.8 


1.8 


9.6 


G29.97-0.05 


9.3 


1.3 (0.7) 


2.6 


4.8 


3.5 


G79.34+0.33 


8.4 


0.7 (0.1) 


2.1 


3.5 


4.9 


G8 1.74+0.59 


16.0 


1.8 (0.2) 


4.7 


3.8 


4.4 



a column density derived from the 850 yum dust continuum flux measured using SCUBA. 

The dust temperature is assumed to be equal to the gas temperature derived from NH 3 (1,1) and (2,2) measurements as given in Table[3] 

b abundance of C ls O relative to H 2 . 

c depletion factor compared to the canonical value of 1.7- 10~ 7 dFrer king et al.l l 19821) . 



Table 6. H CN and HC N line parameters and column density 







H l j CN J = 


1-0, F = 2- 


-1 


HC b N/= 1 


-0 


Source 


f r MB dv a 


Av 


T m { 


^Hl-'CN 


f r MB dv 


Av 






Kkms"' 


km s -1 


K 


10' 2 cm" 2 


Kkms" 1 


kms -1 


10 12 cm" 2 


G10.21-0.31 


0.9 (0.1) 


3.5 


15.4(1.5) 


12.2 (1.1) 


0.2(<0.1) 


1.2 (0.4) 


0.5 (0.1) 


G10.21-0.32 


0.7(0.1) 


2.4 


14.5 (1.7) 


9.3(1.0) 


0.6 (0.2) 


5.6(2.2) 


1.3 (0.6) 


G13. 18+0.06 


3.0(0.1) 


3.7 


17.4(1.3) 


42.7 (2.5) 


0.7 (0.1) 


3.3 (0.4) 


1.9 (0.2) 


G15.01-0.67 


7.3 (0.2) 


3.4 


24.3 (1.5) 


129.5 (7.1) 


2.6 (0.2) 


2.8 (0.3) 


8.6 (0.9) 


G15.01-0.69 


4.8 (0.2) 


2.7 


24.3 (1.3) 


85.2 (4.3) 


1.6 (0.1) 


2.4 (0.2) 


5.3 (0.4) 


G15.03-0.65 


3.6(0.1) 


3.6 


20.4 (2.4) 


55.9 (4.9) 


1.3 (0.1) 


3.6 (0.0) 


4.4(1.0) 




1.6(0.1) 


1.7 


16.6 (2.9) 


21.6 (2.6) 


0.3 (0.1) 


1.7 (0.0) 


0.9 (0.2) 


G23.41-0.23 


1.3 (0.1) 


2.0 


18.8(1.8) 


18.6(1.4) 


0.5 (0.1) 


2.6 (0.5) 


1.3 (0.2) 


G23.42-0.23 


1.4(0.1) 


3.4 


19.0 (2.0) 


20.3 (1.7) 


0.4 (0.1) 


2.6 (0.4) 


1.1 (0.2) 


G28.34+0.06P2 


1.0(0.1) 


2.4 


15.7 (2.0) 


13.8 (1.2) 


0.7 (0.1) 


3.7 (0.5) 


1.8(0.2) 


G35. 19-1.73 


1.9(0.1) 


2.4 


17.0 (0.9) 


27.1 (1.3) 


0.6 (0.1) 


2.0 (0.3) 


1.4 (0.2) 


G79.34+0.33 


0.8 (<0.1) 


1.5 


13.7(1.0) 


10.0 (0.5) 


0.2 (<0.1) 


1.2 (0.2) 


0.5 (0.1) 


G8 1.74+0.59 


3.2 (0.1) 


2.8 


18.4(1.1) 


46.7 (2.3) 


1.1 (0.1) 


2.8 (0.2) 


3.0 (0.2) 



a Integrated line intensity for the H 13 CN main group alone. 
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Table A.2. NH3 line parameters with uncertainties (in brackets) from the hyperfine/Gaussian fits. 







NH 3 (1,1) 






X TT T ^ \ 

NH 3 (2,2) 


Source 


Tmb' 


A h 

Av" 


T c 


Tmb 




(K) 


kms" 1 






(K) 


G8. 13+0.25 


1.9 (0.4) 


1.5 (0.1) 


1.5 


(0.3) 


1.0 (0.2) 




2.4 (0.4) 


1.3 (0.1) 


1.9 


(0.3) 


1.3 (0.2) 


G8.68-0.37 


5.3 (0.5) 


1.1 (0.0) 


5.0 


(0.4) 


2.3 (0.3) 




5.9 (0.4) 


3.0(0.1) 


2.6 


(0.1) 


3.5 (0.2) 


G8.71-0.37 


3.5 (0.4) 


1.3 (0.0) 


4.2 


(0.3) 


1.7 (0.2) 


G10.15-0.34 


1.1 (0.1) 


4.1 (0.3) 


1.1 


(0.5) 


0.7 (0.1) 


G10.21-0.31 


2.9 (0.4) 


1.9 (0.0) 


3.9 


(0.3) 


1.8 (0.1) 


G10.21-0.32 


3.4 (0.5) 


1.7 (0.0) 


3.5 


(0.2) 


1.8 (0.3) 




4.6 (0.4) 


1.9 (0.0) 


4.6 


(0.2) 


2.7 (0.1) 


G10.61-0.33 


2.0 (0.2) 


1.9 (0.1) 


1.8 


(0.3) 


1.1 (0.1) 


G11.11-0.12P1 


5.1 (1.0) 


1.3 (0.1) 


3.5 


(0.4) 


1.9 (0.3) 


G11.11-0.12P3 


3.0 (2.2) 


2.5 (0.1) 


2.4 


(0.2) 


1.0 (0.2) 


G11.11-0.12P4 


4.3 (0.9) 


1.5 (0.1) 


3.1 


(0.6) 


0.9 (0.4) 


G12.19-0.12 


0.4 (0.0) 


3.9 (0.2) 


1.8 


(0.3) 


0.1 (0.0) 


G13.18+0.06 


5.4 (0.5) 


2.4 (0.0) 


3.3 


(0.1) 


3.7 (0.1) 


G 15. 01 -0.67 


5.0 (0.4) 


2.9(0.1) 


0.7 


(0.1) 


3.6 (0.1) 


G15.01-0.69 


3.3 (0.2) 


2.6(0.1) 


0.5 


(0.1) 


2.3 (0.1) 


G15.03-0.65 


2.6 (0.3) 


2.4 (0.2) 


1.2 


(0.3) 


1.6 (0.2) 




2.8 (0.3) 


2.6 (0.1) 


1.3 


(0.3) 


1.2 (0.4) 


G18. 17-0.30 


2.9 (0.3) 


1.6 (0.0) 


3.3 


(0.2) 


1.7 (0.1) 


G18.21-0.34 


3.8 (0.4) 


1.5 (0.0) 


2.8 


(0.2) 


2.0 (0.2) 


G19.30+0.07P1 


5.0 (0.7) 


2.0 (0.1) 


2.1 


(0.3) 


2.7 (0.2) 


G19.30+0.07P2 


3.2 (0.6) 


1.7 (0.1) 


2.9 


(0.5) 


1.4 (0.3) 


G23.41-0.23 


3.0 (0.3) 


2.0 (0.1) 


1.2 


(0.1) 


1.6 (0.2) 


G23.42-0.23 


2.2 (0.3) 


1.9 (0.1) 


1.6 


(0.2) 


1.3 (0.1) 




1.3 (0.4) 


1.5 (0.1) 


0.6 


(0.4) 


0.4 (0.1) 


G23.44-0.18 


3.6 (0.2) 


3.6(0.1) 


1.7 


(0.1) 


2.8 (0.1) 


G27.29+0.15 


1.7 (0.2) 


2.5 (0.1) 


1.6 


(0.2) 


0.9 (0.1) 


G27.31+0.18 


1.0 (0.2) 


2.1 (0.1) 


2.1 


(0.4) 


0.6 (0.1) 


G28.34+0.06P1 


2.7 (0.5) 


2.7 (0.0) 


2.0 


(0.1) 


1.3 (0.1) 


G28.34+0.06P2 


2.8 (0.5) 


2.0(0.1) 


2.7 


(0.3) 


1.5 (0.2) 


G29.97-0.05 


0.8 (0.2) 


2.0 (0.0) 


2.4 


(0.2) 


0.3 (0.0) 


G33.71-0.01 


2.0 (0.4) 


2.4 (0.2) 


2.9 


(0.6) 


1.5 (0.2) 


G34.81-0.28 


4.7 (0.5) 


1.4 (0.0) 


1.6 


(0.2) 


2.2 (0.1) 


G35.19-l.73 


4.7 (0.4) 


1.6 (0.0) 


1.8 


(0.1) 


2.4 (0.1) 


G79.34+0.33 


2.4 (0.3) 


1.2 (0.0) 


1.5 


(0.1) 


0.7 (0.1) 


G81.74+0.59 


6.1 (0.6) 


2.2 (0.0) 


1.3 


(0.1) 


3.2 (0.1) 



11 Main beam brightness temperature from gaussian fit. 
b FWHM from the hyperfine fits. 

c optical depth of the main group of hyperfines from the hyperfine fits. 

Note: Five sources clearly have at least 2 velocity components, hence the fit parameters have been determined for each component. Similar 
notations are used for other molecules elsewhere in this paper. The parameter has been fixed for the fit for those entries with no error indicated. 



Pillai et al. 2007: Deuteration and Depletion in High Mass Clumps, Online Material p 3 
Table A.3. NH2D line parameters with uncertainties (in brackets) from the hyperfine/Gaussian fits. 





NH 2 D 85.93 GHz 


Source 


J r MB dv a 


Av b 


r tot L 




— ^ r , 1 — 

Kkms" 1 


km s~ J 




G8. 13+0.25 


0.5 (<0.1) 


1.3 (0.3 ) 


- 




0.7 (<0.1) 


1.8 (0.9 ) 


- 


G8.68-0.37 


- 


1.2 (0.1 ) 


2.2 ( 0.7 ) 


G8.71-0.37 


- 


1.1 (<0.1) 


3.7(0.6) 


G10.21-0.31 


- 


1.2 (0.1 ) 


3.4 ( 1.4) 


G10.21-0.32 


- 


1.4 (0.1 ) 


2.3 ( 1.0) 


G10.61-0.33 


0.3 (<0.1) 


1.5 (0.2 ) 


- 


G11.11-0.12P1 


- 


1.0 (0.1 ) 


3.6 ( 1.5 ) 


G11.11-0.12P4 


0.4(0.1 ) 


1.3 (0.2 ) 


- 


G12.19-0.12 


1.5 (0.1 ) 


3.6 (0.5 ) 


- 


G13. 18+0.06 


- 


2.2 (0.1 ) 


2.0 ( 0.5 ) 


G18. 17-0.30 


- 


1.5 (0.1 ) 


4.4 ( 1.1 ) 


G18.21-0.34 


- 


1.3 (0.1 ) 


2.0 ( 0.8 ) 


G23.41-0.23 




1.7 (0.2 ) 


2.1 ( 1.3 ) 


G23.44-0.l8 


1.4 (0.1 ) 


3.1 (0.2) 




G27.29+0.15 


0.4 (0.1 ) 


1.5 (0.4) 




G28.34+0.06P2 


1.2 (0.1 ) 


1.8 (0.2 ) 




G29.97-0.05 




1.4 (0.1 ) 


3.0 ( 0.6 ) 


G33.7l-0.0l 


0.4 (<0.1) 


1.6(0.3) 




G34.8l-0.28 


0.8 (0.1 ) 


1.7(0.1 ) 




G35.l9-l.73 




1.4 (<0.1) 


2.7 ( 0.4 ) 


G79.34+0.33 




1.0 (0.9) 


1.1(0.1) 


G81.74+0.59 




2.8(0.1 ) 


1.5 (0.3 ) 



a Integrated intensity of the line from the Gaussian fits. 

b FWHM from the hyperfine fits. 

c Total optical depth from the hyperfine fits. 

Note: Similar notations are used for other molecules elsewhere in this paper. The parameter has been fixed for the fit for those entries with no error 
indicated. 
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Fig. A.l. Lower spectrum in each panel: Effelsberg 100m observation of the NH3 (1,1) emission. Upper spectrum in each panel: 
30m observation of the NH 2 D at 85.9 GHz. The NH 2 D is scaled by a factor 5.0 to amplify the emission relative to NH 3 in the 
absolute (brightness temperature) units. 
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LO O lO 
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Fig.A.2. lower spectrum:30m C 18 (J = 1-0). upper spectrum: 30m C 18 (J = 2-1). The velocity range is ±15 kms" 1 of the 
systemic velocity. 
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Table A.4. C ls O (J = 1-0) and (J = 2-1) line parameters 





L. U (J = i— U) 


C <J ( J = L— 1 ) 


Source 


fT MB dv 


[TmbcIv 




Kkms-' 


i r , 1 

Kkms-' 


G8.13+0.25 


15.7 (0.1) 


27.6 (0.5) 


G8.68-0.37 


24.9 (0.9) 


31.9 (1.7) 


G8.71-0.37 


13.5 (0.4) 


5.9 (0.7) 


G10. 15-0.34 


15.6 (0.3) 


37.4 (0.5) 


G10.21-0.31 


14.9 (0.2) 


19.0 (0.4) 


G10.21-0.32 


27.5 (0.3) 


33.2 (0.6) 


G10.61-0.33 


10.6 (0.3) 


14.7 (0.5) 


G11.11-0.12P1 


6.1 (0.1) 


5.5 (0.4) 


G11.11-0.12P3 


5.1 (0.3) 


5.1 (0.5) 


G11.11-0.12P4 


5.1 (0.2) 


4.8 (0.5) 


G12.19-0.12 


12.0 (0.2) 


20.8 (0.4) 


G13. 18+0.06 


9.7 (0.5) 


20.9 (0.8) 


G15.01-0.67 


31.0 (0.2) 


91.7 (1.0) 


G15.01-0.69 


16.2 (0.2) 


53.7 (0.5) 


G15.03-0.65 


19.8 (0.2) 


65.5 (0.6) 


G18. 17-0.30 


13.1 (0.7) 


26.5 (1.3) 


G18.21-0.34 


7.3 (0.4) 


12.0(1.3) 


G19.30+0.07Pl a 


1.1 (0.4) 


- 


G19.30+0.07P2 


4.1 (0.4) 


6.6 (0.7) 


G23.41-0.23 


35.2 (0.8) 


44.5 (1.0) 


G23.42-0.23 


31.5 (0.3) 


43.1 (0.5) 


G23.44-0.18 


31.1 (0.5) 


43.4 (0.6) 


G27.29+0.15 


12.4 (0.2) 


16.6 (0.5) 


G27.31+0.18 


4.9 (0.2) 


6.1 (0.5) 


G28.34+0.06P1 


11.3 (0.4) 


6.3 (1.5) 


G28.34+0.06P2 


14.4 (0.5) 


19.1 (0.7) 


G29.97-0.05 


12.9 (0.7) 


18.2 (0.9) 


G33.71-0.01 


7.4 (0.2) 


9.0(1.3) 


G34.81-0.28 


5.4 (0.3) 


6.6 (0.6) 


G35.19-l.73 


8.4 (0.2) 


15.8 (0.8) 


G79.34+0.33 


7.8 (0.1) 


9.9 (0.2) 


G8 1.74+0.59 


14.9 (0.2) 


32.4 (0.4) 



j r MB dv is the integrated intensity over ±5 km s 'of the brightest C 18 (J = 1-0) peak. 
a C ls O (7 = 2-1) data suffered from baseline ripples. 
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Fig. A.3. lower spectrum: 30m observation of H 13 CN . upper spectrum: 30m observation of HC 15 N. The HC 15 N is scaled by a factor 
2 to amplify the emission relative to H 13 CN in the absolute units. The velocity range is ±25 kms~' of the systemic velocity. 



